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Abstract. We study a class of Hamiltonian deformations of the massless Einstein- 
O ! Klein-Gordon system in spherical symmetry for which the Dirac constraint algebra 

closes. The system may be regarded as providing effective equations for quantum 



. gravitational collapse. Guided by the observation that scalar field fluxes do not follow 

metric null directions due to the deformation, we find that the equations take a simple 



form in characteristic coordinates. We analyse these equations by a unique combination 
of numerical methods and find that Choptuik's mass scaling law is modified by a mass 
gap as well as jagged oscillations. Furthermore, the results are universal with respect 
OO . to different initial data profiles and robust under changes of the deformation. 
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1. Introduction 

Einstein's theory of general relativity (GR) predicts the existence of black holes. 
According to the Penrose- Hawking singularity theorems they can form when matter 
undergoes gravitational collapse. However, the singularities inside black holes are 
unwelcome from a physical point of view as they render spacetime incomplete. We 
have no knowledge about what happens near a singularity. This is a regime where 
quantum gravity effects are expected to come into play. 

The canonical formulation of GR is one approach to formulating a theory of 
quantum gravity. It can be traced back to the work of Arnowitt, Deser, and Misner 
(ADM) [lj who gave the first Hamiltonian formulation. This was used by Wheeler [2] 
and DeWitt |3j to formulate the quantization program. In practice, there are many 
technical and conceptual difficulties in this approach [3], and so far no complete theory 
is available. 

In the absence of a full theory it is important to ask whether there are relevant 
physical situations with symmetries where the canonical quantization can be completed. 
Technically the easiest is that of a homogeneous cosmology, in which case GR reduces 
to a system with finitely many degrees of freedom. Such reductions have been studied 
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since the 1970s using the Wheeler-DeWitt approach [5], and again in the recent past 
using loop quantum gravity (LQG) methods [6J. 

The next level of problem, beyond quantum mechanics, is dimensional reduction 
to a tractable field theory. Two important problems of physical interest that fall into 
this category are that of a cosmology with inhomogeneities and that of the spherically 
symmetric gravitational collapse. Our interest in this paper is in the second problem. 
We address the question of how to formulate deformations of Einstein's equations in 
spherical symmetry and discuss how they can be interpreted as providing "effective" 
quantum gravity corrections using some input from LQG. 

Classical spherically symmetric collapse equations have been extensively researched. 
Using Schwarzschild coordinates [TJ, Christodoulou showed [8] that there are two classes 
of initial data for the collapse. For one of them ("weak data") matter bounces at the 
centre of the coordinate system and evolves back towards radial infinity. The other class 
of data, however, gives rise to a black hole. This was numerically confirmed not only by 
Choptuik [7J but also by Goldwirth and Piran [9], who used infalling null coordinates 
|10] . For a comparison of the two approaches see [11]. The next important development 
in the history of classical spherically symmetric gravitational collapse was by Choptuik 
[12] . Among other things, he found numerically a scaling law which relates the amplitude 
of the scalar field to the mass of the black hole. This was later reproduced by Garfinkle 
[T3] based on collapse equations in double null coordinates. Such a form of the equations 
was also used in [T3] to find a similar scaling law with a nonzero cosmological constant 
and in higher dimensions. 

The investigation of quantum corrected spherically symmetric collapse was the 
natural next step to take. In |15j it was shown that certain corrections motivated 
by the LQG program introduce a mass gap in the mass scaling law. This indicates 
that zero mass black holes are no longer contained in the solution space. Using similar 
corrections in Painleve-Gullstrand coordinates |16| . Ziprick and Kunstatter [T7 ] fT8 ] also 
found the mass gap, and importantly, were able to observe the evolution beyond the 
horizon. 

In none of the works mentioned in the previous paragraph have the modifications 
been directly introduced in the context of a Hamiltonian formulation. Rather, the 
corrections were inserted into the evolution equations or a covariant action (in [18j this 
is the so-called dilaton action for GR). Also, at least in [T5l [T7] . it has not been checked 
whether the modified equations are consistent with GR's symmetry group, that is the 
diffeomorphisms. It is therefore interesting to ask if it is possible to arrive at effective 
equations representing quantum gravity corrections from a Hamiltonian formulation and 
that are such that the Dirac algebra closes. This forms part of the motivation for the 
present project. 

Following work on loop quantum cosmology by Bojowald and his collaborators |19] . 
where a class of anomaly-free effective constraints are presented, Reyes |20j obtained a 
closed constraint algebra in spherical symmetry using the connection-triad variables 
of LQG. However, the resulting equations are sufficiently complicated as to make a 
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numerical investigation quite difficult. 

The equations we derive and analyse in this paper follow from the same general 
ideas but are obtained directly in the ADM formulation. The main obstacle faced with 
the connection-triad variables is overcome to yield a remarkably simple form of the 
equations that can be compared directly to the classical double null formulation. Along 
the way we note a useful physical insight crucial for the form of the equations we obtain: 
the null coordinates of the metric do not coincide with the characteristic lines of the 
scalar field due to the deformation. This is what may be expected to be necessary for 
singularity avoidance, since it signals a violation of the dominant energy condition. 

The paper is structured as follows. In the next section we present the canonical 
theory and introduce a class of deformations that are to model quantum corrections. 
In Section [3] we derive the corresponding canonical equations of motion for the 
Schwarzschild gauge. The transformation to characteristic coordinates, which yields 
the simple equations, is given in Section HI In Section [5] we choose a specific function 
that deforms GR. The numerical scheme to solve the resulting equations is defined in 
Section [6l The definition of black holes requires some commentary, which is provided in 
Section [71 We proceed in Section [8] by posing the questions that we want to answer in 
the remainder of the paper. Finally, Section [9] contains the results and Section [TUl some 
concluding remarks. 



2. Canonical formulation and its deformation 

In metric ADM variables the line element for a spherically symmetric spacetime can be 
defined by 

ds 2 = — [N 2 — {N r A) 2 ]dt 2 + 2N r A 2 dtdr + A 2 dr 2 + B 2 dQ 2 , (1) 

where N is the lapse function and N r the radial component of the shift vector. If, to 
this metric, we minimally couple a massless scalar field with matter Lagrangian density 
— a/— \g\ g^d^dycf)/ (8tt), we get the canonical actional 

ft\ poo 

S= / (paA + pbB + p^-H) dr dt. (2) 
Jto Jo 

The Poisson brackets are 

{A(t,r),p A (t, s)} = {B{t,r),p B {t,s)} = {<f){t,r),p^t, s)} = 5{r, s) 

and the total Hamiltonian density H = NH± + N r, H\\ r consists of the Hamiltonian 
constraint densitjj|| 

H ± = 7^(paA - 2 PB B) - ^[A'B 2 ' - A(2BB" + B' 2 ) + A 3 ] 



ir B 



(3) 



2AB 2 2A 

fThe boundary terms can be ignored here. 
\ We set G N = 1 and c = 1 . 
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and the diffeomorphism constraint density 
H\\ r = -p' A A + p B B' + prf'. 

From this well known setup [21] we can introduce a classical deformation that is 
designed to describe a certain type of quantum effect. As noted above, the way we 
arrive at the correction is based on the ideas of Bojowald and his collaborators [19j, 
which were also applied by Reyes [20] ■ One of the differences here is that we are not 
using the connection-triad variables [22] . 

The insight motivating the deformation is that the inverse powers of A and B in 
(J3J) can be conveniently quantized using the Thiemann prescription for the inverse triad 
operator (23jJtl Doing so, we can turn H± = / °° NT-L± dr into an operator H±. If we 
proceed by computing the expectation value of H± for certain classes of statedjj, we get 
an effective H e ± . The afore cited work by Bojowald and others suggests that it can be 
of the form Hj_ = J °° iV"H]_ dr, where the effective Hamiltonian constraint density is 

3 

WL^QqH® (4) 



with 



n ( l ] = ^- 2 ( PA A - 2p B B) - ^~ 2 [A'B 2 ' - A{2BB" + B' 2 ) + A 3 

nj{2) P± n/(3) _ B 2 (j)' 2 

ft I = - , „ n , ft 



2AB 2 ' ± 2 A 

In this definition, the Qu\ depend only on B and the describe, in increasing 

order of i, the extrinsic curvature plus the Ricci scalar of constant parameter time t 
hypersurfaces, the kinetic energy of the scalar field, and its gradient energy. The form 
of 01]) implies that the deformation is entirely contained in the functions Qu). Moreover, 
by restricting the Qu\ to depend only on B, modifications due to terms involving inverse 
powers of A are disregarded. Lastly, the term stemming from the extrinsic curvature and 
the one from the Ricci scalar are treated equally even though their overall dependence 
on B is different. 

In spite of these concerns, there is a very good reason to investigate the physical 
consequences of (j3J). Namely, if we impose the condition 

Q\x) = Q(2)Q(3) (5) 

it can be shown (as in [20]) that the Dirac constraint algebra closes. There are no 
anomalies if fl5]) is true. This has far reaching consequences. Among them is that all 
degrees of freedom remain to be captured by the given canonical variables. Furthermore, 
it implies that energy is conserved, which is attractive from a physical point of view. 

f We note that this operator is not the only source for quantum corrections; the other is the curvature 
written using holonomy operators. It is not known how to get a deformed, anomaly-free effective system 
that includes this correction in general, but some special cases have been studied [201 124) . 

|See for the case of cosmological models. 
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3. Equations in Schwarzschild coordinates 

To get an idea of the effect of the quantum corrections introduced in the previous section, 
we implement the second class gauges Q» r = r — B = and Q± = p^ = in the given 
order to freeze out some of the degrees of freedom [26J. Concretely, the consistency 
condition Q\\ T = gives N r = Npa/t and the equation %\\ r = yields 

Pb = Pa A - P</>4>'- (6) 

Note that the gauge choice Q\\ r = renders the nondynamical by turning them into 
functions that solely depend on r. Proceeding with the gauge fixing, the choice of Q± 
forces iV r to vanish, which not only turns the line element ([TJ into the familiar form 

ds * = _AT2 dt 2 + A 2 dr 2 + r 2 dn 2 (?) 

but also, from fl6]) and the definition of ps, implies the trivial equation 

A = Q { i)Np^^ 
r 

The consistency condition Q± = leads to 
N' _ A' A 2 - 1 Q' {1) 

N A + r Q {1) [y) 
and solving T-L\ = for A gives 

A' I- A 2 1 



A 2r ■ SS^Wm^ + ^O. (10) 
Finally, Hamilton's form of the scalar field equation is given by 

If we set all Qu) equal to 1, the deformation is switched off. This allows us to 
recover in (|S1)- (jlip the familiar collapse equations in Schwarzschild coordinates. This is 
the form of the equations that has been heavily investigated analytically and numerically 
(see [27J for an extensive list of references). 

However, to reproduce Choptuik's numerically found results in these coordinates 
is a rather complex task. The conceptual reason is that near r = 0, where black holes 
are expected to form, a numerical code has to allow for a very high resolution but 
as we move away from the origin there is no need for the same accuracy. Trying to 
keep the computing time low, Choptuik decided to use the adaptive mesh refinement 
algorithm of Berger and Oliger |28| . Unfortunately, and despite its beauty, the method 
is rather involved so that the conceptual problem is effectively replaced by a practical 
one. Because of this we seek a different form of the deformed equations, one that is 
more user-friendly. 
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4. Equations in characteristic coordinates 

Christodoulou [10J showed that Einstein's field equations for the model defined by (j2J) 
can be given in a very compact form if infalling null coordinates (u, r) are used. Some 
years later, Garfinkle [13] derived the same equations in double null coordinates (u,v). 
Namely, in terms of a line element of the form 

ds 2 = -W 2 dudv + r 2 dQ 2 , 

Einstein's field equations are given by |29] 

W 2 = -2d u d v r 2 , 

2d u r^- = d 2 u r + r(<9 u 0) 2 , 

dW W 

2d v r^- = d 2 r + r(^0) 2 , 
d u rd v (p + d v rd u (f> = -rd u d v (fi, 

and due to the Bianchi identities the second and the third equation are redundant [10J. 
Dropping the second, Garfinkle used the reparametrization W 2 = 2d v rF to arrive at 
the evolution equations 

d ur = ~, <9„$ = ^(/-F)(0-$), (13) 

where 

r 

^F = ^F(0-$) 2 , (14) 
r 

r 

These equations are of the same compact form as those presented by Christodoulou for 
infalling null coordinates. 

Equations (Tl3|) and (I14p have several advantages compared to those in Schwarzschild 
coordinates ()8l)-( lTT|) (with the deformation switched off). The most important is that 
in double null coordinates (u, v) we have d u r < 0, which follows from the first equation 
in ( fT2l [29J. This implies that the computational grid becomes smaller as r = is 
approached so that a higher accuracy can be naturally obtained by inserting additional 
points into constant u slices [TBI [HI [15]. An adaptive mesh refinement algorithm in the 
spirit of that used by Choptuik is no longer necessary. We therefore intend to transform 
our deformed equations to double null coordinates. However, as we will see, this is not 
enough to obtain simple equations due to a null characteristic mismatch. 

To simplify the following expressions, we introduce the functions Q = (Q(2)/Q(i)) 1 ^ 2 
and s ee A/(Qn\N) as well as the matter variables 

x = ( s 0-0')/Q, ^ee(4 + 0')/Q. 
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This allows us to write the trivially solved Equation (|S]) in the form 

^ = -^(x 2 -V> 2 ) (15) 

and the remaining Equations f[9"j)- (lTTj) as 

s' I -A 2 A' I -A 2 



s r A 2r 4 

sx + x' = ~ Hr/s)'x + ln(r/Q)V, (16) 
sip — ip' = — ln(r/Q)'x + ln(r/s)'ip. 
The last two equations imply that we can denote by c±t? M = sd t ± d r the characteristic 



directions |30j of the scalar field equation. By definition, the variable x describes the 
change of (p along c_ with speed — s, and ip does so along c + with speed s. Therefore, 
the last two equations in (1T6|) tell us how a consecutive change of along c_ and c+, 
and vice versa, looks like. 

Note, however, that the characteristics do not correspond to the null directions 
Q(i)sd t ± d r of the metric ((7j). In fact, we have 

<W'-±-'-±)' J = -^^ La ^ (1') 

V(i) 



which vanishes if and only if Qm = 1. This means that the fluxes of the massless scalar 
field do not follow null directions so that the theory under investigation is not GR. 

The important implication of (TTT)) is that only if Q 2 ^ < 1 are we dealing with 
matter of which the fluxes are timelike or null. If Q%\ > 1 the scalar field propagates 
along spacelike directions, giving a violation of the dominant energy condition (see also 
[18J). In this case we can hypothesise that the formation of singularities can be avoided. 

What kind of matter are we dealing with? In this context, it is interesting to note 
that there is no scalar field potential that can be added to the Lagrangian matter density 
to reproduce (TI5]) and f|T6l) by means of Einstein's equations. This follows from the fact 
that a potential not only gives a contribution to the scalar field equation but also to the 
consistency condition and the constraint equation. However, in the given variables the 
latter two are effectively unaffected by the quantum corrections (Q only appears in the 
last two equations in (TL6]) ). 

The equations for x an d ip in ( fT6l) suggest that we introduce coordinates u and v 
adapted to the characteristics c±. We therefore set 

D u = -sd t + d r , D v = sd t + d r , (18) 

where D u = {d u r)~ l d u and D v = (c^r)"" 1 ^. This gives 

X = -D u <p/Q, V = D v <f>/Q. 

That is, since d u r < and d v r > (see below), the coordinates u and v respectively 
parametrize future pointing infalling and outgoing characteristics if Q 2 ^ < 1. For (fT8|) 
we chose d u t = —sd u r and d v t = sd v r, which implies the coordinate relation 

- W 2 dudv = -Q 2 {1) N 2 dt 2 + A 2 dr 2 (19) 
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if the function W is parametrized by 
W 2 



... , =A*. (20) 

Since the characteristics c± to which we adapted the coordinates u and v do not agree 
with the null lines of the physical relevant line element we of course cannot expect 
u and v to correspond to the double null coordinates of this metric. As is shown by 
ffT9]) . u and v correspond to the double null coordinates of that metric, for which the 
null directions agree with the characteristics. 

If we factor out A 2 in (1T9]) . we can show from (120]) that s'/s = —d u d v r/(d u rd v r). 
This, together with the first equation in ( fT6l) . implies 

W 2 = -2d u d v r 2 . (21) 

A comparison with (|12|) shows that the deformation has no impact on this equation. 
From ( 121"]) it follows that we again have d u r < [29] and since ( 120]) implies that 
d u rd v r < we get d v r > 0, as desired. Constructing linear combinations sA ± A' 
from ( fl~5]) and ( 1T6]) . Equation ( 120]) gives 

2<V^ = ^ + r(d u 0/Q) 2 , 

(22) 

2<%r^- = ^r + r(^/Q) 2 . 

Since we are no longer dealing with Einstein's field equations, the Bianchi identities 
cannot be used to show a redundancy in f ]2"2"]) . Luckily, the trivially solved Equation 
(1T5]) can be used here. Note that the equations in (122]) do not differ very much from the 
corresponding ones in ( TT2"]) . Finally, the equations for \ an d 4> i n (116]) give the scalar 
field equation 

d u (r/Q)d v <j ) + d v (r/Q)d u <f ) = -(r/Q)d u d v <j) (23) 

and it is this equation, in which the quantum corrections are most prominent Jjj 

As expected, the speed s does not appear in any of ( ]2T]) - (l23]) . It cannot be "seen" 
along the characteristics. Therefore, the only quantum corrections left are encoded in Q. 
This implies that we can switch them off not only by setting Q(i), Q(2), and therefore 
Q(3) equal to 1 but also by means of the relation Q(i) = Q(2) — Q(3)- This is not 
very surprising since, as we have seen, the gauge Qj_ turns the into nondynamical 
functions and if they are the same, (J2]) and @ show that we can effectively remove them 
by reparametrizing the lapse function. The class of quantum corrections where all 
agree are therefore uninteresting from a physical point of view. 

To arrive at a form of the evolution equations that resembles the one given in (1T3]) 
and (JHJ) we set W 2 = 2d v rF. Dropping the first equation in ( 122]) . we can focus on the 
consistency and the constraint equation 

d v rF = -d u d v r 2 , d v r^f = r(^/Q) 2 , (24) 

fit is also apparent that one could obtain a modified dispersion relation from this equation, which 
is a topic of much recent interest. 
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respectively, as well as on ( I23p . The key now is to realize that the relevant radial 
coordinate in (|24|) is r, whereas in (123]) it is r/Q. To remove this asymmetry we define 

which can be seen as the result of exchanging Q/ 2 ) with r 2 Q( 2 )- The function g 
behaves like Q in the sense that it is nondynamical with respect to the Schwarzschild 
coordinates (t,r). In terms of the characteristic coordinates (u, v ) its dynamic is 
implicitly determined by that of r. Finally, since q = implies D u q = D v q = q', 
we can arrive at the sought after evolution equations. They are 



d u r = — , <9 U $ 

2 2r 



- *), (26) 



where 



^ F= a^ (27) 

r 
r 

A comparison with (TT3j) and (TT4j) shows that the deformation is elegantly captured by 
relatively minimal modifications. 

It is now straightforward to produce an evolution scheme for these equations. This 
requires only one additional input, which is the selection of the function q such that 
it gives modifications to the classical behavior in regions of sufficiently large energy 
density. 



5. Modifications 

The definition of 7-L\ in (J3J) suggests that the functions Q^(B) modify the factors 1/B 2 . 
Namely, according to |31[ [T5] . such modifications can be explicitly constructed from 
the expectation value of the operator 1/B on the so-called polymer Hilbert space [32J. 
Originally, this method was motivated by the use of the polymer quantization in the 
context of LQG and loop quantum cosmology, and the result that we want to take 
serious is that it can remove the functional singularities of 1/B 2 . That is, 1/B 2 can be 
"regularized" on the real line and we will choose the to do precisely this. In this sense 
we can interpret these functions as describing effective quantum gravity corrections. 

Since the gradient energy is not inversely proportional to B 2 , we can set 

Q(3) = 1. The condition ([5]) yields then Q 2 ^ = Q( 2 ) and we are left with one independent 
Qu). If we now go back to the gauge fixed sector where B = r, the function q from ( |25l) 
becomes 
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Figure 1. The function q = yl + r 2 ' modifies the GR curve q = r. Their difference 
vanishes in the limit r — > oo. 



Because of f lT7|) this means that the characteristics u, v £ R of the scalar field equation 
are timelike wherever g > r and spacelike for q < r. 

Given our intention to think of Q(i) as a means to regularize 1/B 2 = 1/r 2 , we require 
1/q < oo for r — >■ 0. We also demand that the impact of g vanishes in the asymptotically 
flat region, where energy densities and, thus, modifications due to effective quantum 
gravity corrections should be negligible. The remaining conditions on q are that it is 
strictly increasing and positive as, otherwise, either the term q"/q' in ( 1261) or the term 
q'/q in (j27|) lead to singular evolution equations. A valid choice is therefore 

q = V\ 2 + r 2 \ (28) 

where we can use the length scale A to give a more precise definition of the condition on 
q in the asymptotic region. Namely, we require |r* — gj <C A for r A. Note that (1281) 
implies that the characteristics u, v £ K. are everywhere timelike. 

Having specified q, we can introduce dimensionless variables u, v, r, and q by 

{u, v, r, q} i — > {Au, Av, Ar, Aq}. 

The equations of motion (126]) and ( )27|) are now dimensionless. The function q is now 
(see Figured]) 

q = Vl + r 2 ' (29) 
and we can move on to the numerics. 

6. Numerical scheme 

The boundary conditions are r(u, u) = [29] and that $ can be expanded in a series in 
v in a neighborhood of v = u. That is, we assume the existence of an expansion 

$(u, v) = $(u, u) + (<9 v $)(u, u)(v - u) + 0((v - u) 2 ), (30) 

which was also done in [13]. The boundary conditions imply 

F(u, u) = 2(cV)(u, u) = /(u, u), («9 v F)(u, u) = = (d v f)(u, u), 
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for the metric function F and the auxiliary field /. The first relation is a consequence 
of the field equation d v rF = — <9 u <9 v r 2 together with (<9 u r)(u, u) = — (<9 v r)(u, u). The 
second relation is only given as it is used in the implementation of the v-integrator 
specified below. Regarding the boundary conditions for the matter fields, we have to 
have 0(u, u) = $(u, u) in order to get evolution equations that are regular at r = 0. 
Finally, we need to know (<9 v 0)(u,u) and (<9 u $)(u,u) as they enter the setup of the 
v-integrator and the u-integrator, respectively (see below). They both depend on the 
choice of q and for ( 129]) we have 

(d v 0)(u, u) = 0, («9 u $)(u, u) = (<9 v $)(u, u). 

The initial data for r is specified by 

r(0,v) = ^ (31) 

which implies F(0,0) = 1 [29j. While the boundary condition r(u, u) = removes the 
possibility of a reparametrization u i— > u(u), ( 13TT) makes v i— > v(v) impossible. The 
remaining initial data is for the matter sector and is specified by a choice of 0. Guided 
by [12] and others, we consider the profiles 

0(0, v) = O -^ exp (Jl^l) , (32) 

0(0, v) = 0o tanh ■ (33) 

The dependence on v in (132]) differs from that in the mentioned literature, where it is 
either o v 2 exp(. . .) or o v 3 exp(. . .). The factor v 3 /(l + v 3 ) allows us to choose v 3> 1 
and, thus, to set up the initial data in the regime where the impact of q is minimal. 
That is, we do not have to worry about a black hole being already present on the slice 
u = 0. However, since v 3 / (1 + v 3 ) = l + 0(v -3 ) for v — > oo, it seems likely that this factor 
disallows us to find a critical Vq 3> 1 for which Choptuik's mass scaling law applies (see 
Section E]). 

Our interest in the remainder of this paper is mainly in the initial data profile 
given in ( 132]) . The profile in (133]) is implemented in only one simulation and for the 
completeness of the analysis. It is important to note that (132]) is consistent with the 
boundary condition 0(u, u) = $(u, u) since it implies (<9 2 0)(O,O) = 0. The latter is 
necessary because we have to solve $ = + 2q<9 v 0/q' for $ on the initial slice, and 
the q in (129]) is such that q'(O) = 0. The initial data in (133]) satisfies the condition 
(<9 2 0)(O,O) = the better the lar ger Vo is. Finally, given the form of 0(0, v) in ( 132]) 
and the above comments on v , the parameters of interest are 0o and w, that is the 
amplitude and the width of the Gaussian-like wave packet, respectively. 

To solve ( 126]) numerically, we need a u-integrator. Since we are effectively dealing 
with a modified wave equation, we choose the Richtmyer two-step Lax-Wendroff method 
|33j. It is particularly well adapted to the computational domain, of which an illustration 
is given in Figure [2j Note that the domain becomes smaller as u increases. Garfinkle |13| 
made use of this by injecting additional points into slices of constant coordinate time 
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u 



t 



V 



U = V 




V = V, 



u = 



r 



max 



Figure 2. The numerical domain of integration. For q ^ r the lines u, v 
coincide with the null directions of the metric. 



S R do not 



u, whenever one half of all the points disappeared beyond the boundary at r = 0. This 
is a simple method to achieve an increase in the resolution as the evolution progresses. 
However, we do not use this method and derive all results on the basis of a uni-grid code. 
Finally, we choose the well known second-order Runge-Kutta method as the v-integrator 
to solve ( 1271) . 

The resulting code is second-order accurate [7J EI] in the variables r, $, and but 
only first-order accurate in F and / |35j. An explanation for the lower accuracy of F 
and / is that they are derived not only from r and $ but also from <f), which itself is 
derived from r and $. That is, the inherited error can be expected to be much bigger 
for F and /. According to [36], an increased inherited error is a common feature of 
characteristic codes anyway. However, since our code reproduces Choptuik's classical 
findings, the lower accuracy of F and / is not a problem. 

7. Black holes 

Since we are interested in the impact of various functions q on the black hole formation, 
we need to define a criterion that specifies black holes in characteristic coordinates 
(u, v). With respect to the line element —N 2 dt 2 + A 2 dx 2 + r 2 dVt 2 such a criterion can be 
formulated with the help of the null expansion 



of a sphere of radius r embedded in a slice of constant coordinate time t. Namely, a local 
minimum #™ m that is such that 9™ m — > as t increases, indicates the presence of an 
apparent horizon and, thus, the presence of a black hole. Note that in the coordinates 
(t, r) we have no access to the expansion 9_ of a congruence of infalling null geodesies 

Emm Eg fl 

f This is also true for double null and characteristic coordinates. 



e. 




+ = 
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Since 8 + is a scalar, the corresponding expression in characteristic coordinates (u, v) 
is simply 



A good practical criterion for black hole formation in our model of deformed GR would 
therefore be to look for coordinates (uah ; v ah) f° r which #+(uah ; v ah) is smaller than 
some user-defined threshold 6+ G R+. The code could then be stopped when 9 + < + 
and the mass Mah = ^ah/2 = r(uAH> v AH)/2 of the apparent horizon computed|j] 
However, since in characteristic coordinates (and in double null coordinates as well) the 
code terminates "by itself" when an apparent horizon forms (the coordinates become 
singular) there is no need for the criterion 9 + < 0+. We simply keep a record of the 
values of #™ m and use the last one before the termination to get Mah- 

For q = r the characteristics follow the null geodesies and the mass Mah is identical 
to the black hole mass Mbh- However, at points where q < r the slices u G R are spacelike 
and they are timelike where q > r. It is then not clear from first principles whether we 
can set Mah = Mbh- It is interesting that for q ^ r the apparent horizons tend to form in 
the asymptotic region v ^> 1 of the computational domain but, unfortunately, this does 
not mean that they can be expected to coincide with the event horizon. The reason for 
this is that r remains small there, that is r < 1, which is exactly where the modifications 
are active. To address this issue, that is to understand the relation between Mah and 
Mbh for q ^ r, we should really solve the geodesic equation over the entire spacetime 
(see |39j in this context). On the one hand, it is unfortunate that this cannot be done 
in the given coordinates. On the other hand, we do not have to be too concerned about 
the causality implied by space- and timelike slices as the results we get for Mah are very 
robust in that they appear insensitive to the choice of q (see the various plots below 
and, in particular, Section [9.1. 51) . 

8. Questions of interest 

The questions we are interested in refer to Choptuik's mass scaling law [12| [27] . As far 
as the parameter O is concerned, it reads 



where 0q is a critical amplitude in the sense that initial data with O > 0o leads to 
a black hole, whereas smaller values do not. The value of the exponent 7 in (134j) is 
approximately 0.374 [ID] and Choptuik provided evidence for it to be universal. That 
is, the same value can be found for all members of a family of initial data profiles such as 
the members ( 1321) and ( 1331) of the family of massless scalar field profiles. Furthermore, 

fin (t, r)-coordinates one sets A 2 = (1 — 2M/r)~ 1 and black hole formation is signaled by A — > 00, 
that is M -> r/2. 




Mah oc (0 o - 



(34) 
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Choptuik showed that for (I32p (with v 3 / (1 + v 3 ) replaced by v 3 ) there is a critical value 
w* for which (134"|) reads 

Mah oc (w* - w) 7 . (35) 

The exponent 7 can once again be approximated by 0.374. 

The form of ( ]3~4"]) and (135]) implies that there is a naked singularity with Tah — 
corresponding to the parameters 0q an d w*. It can therefore be expected that quantum 
gravity corrections change the relations f l3~4"|) and (135]) . Indeed, this expectation was given 
more credibility by the work in [T5l H7] . where it was shown that similar corrections to 
those we consider lead to a mass gap. This means that a finite mass is required to form 
a black hole, which suggests that the naked singularity can no longer be found. 

The questions we are interested in are as follows. 

1) Do (134]) and (135]) remain valid in the presence of various functions q? 

2) Do the profiles given in (132]) and (133]) lead to different results? 

3) Is there a mass gap? 

Another interesting point concerns the oscillations that should be incorporated in 
the mass scaling laws ( 13~4"]) and ( 135]) |4T] HQ]. The following question refers to these 
oscillations. 

4) What happens to the fine structure of ( |34j) and ( 135]) in the presence of various 
functions q? 

We will first address the questions 1), 2), and 3), and then turn our attention to 
question 4). 

9. Results 

The figures in the remainder of this paper are the result of simulations performed at a 
resolution A = 2~ 8 ~ 4 • 10 -3 . Furthermore, in all simulations but the ones leading to 
the GR points, we choose a grid with v max = 400 ± 100 and an initial position vo = 80, 
which on the initial slice is significantly larger than the extent of the modified domain 
v < 2. In the GR runs a smaller grid with v max = 240 is used but because of our intent 
to make a comparison with the modified GR runs, we do not alter Vo- 

With the exception of the computations needed to analyse the validity of (135]) . 
where the width w is varied and the amplitude 0o is fixed, we set w = 32 and varied 0o- 
The implemented initial scalar field profile is always that given in (132]) . except for the 
simulations addressing the question 2) above, where it is the one in (l33j) . Finally, in all 
cases the point on the very left-hand side agrees with the critical parameter (0q or w*) 
by 14 decimal places. For all of the following mass scaling plots the simulations have 
been executed on the Atlantic Computational Excellence Network. 
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Figure 3. An illustration of the impact of q = %/ 1 + r 2 ' on Choptuik's mass scaling 
law. Notice the mass gap, the oscillations, the overlap for log 10 ((/>o — <f>o] ^ — 4, and 
the scaling with respect to A. The slope 7 ~ 0.376 of the straight line about which the 
GR data oscillates was determined by a least squares fit over the interval [—13, —3]. 

9.1. The mass scaling law 

Figure [3] shows that in the presence of q = y/1 + r 2 ' Choptuik's mass scaling law is 
significantly modified. The most obvious features of the new data points are a mass gap 
M g , an oscillatory behavior for log lo (0o — < —4, an overlap with the GR points for 
log lo (0o ~ 0o ) ~ — 4, and a scaling with A (note that Mah is denoting Mah/A). Finally, 
we will argue further down in favour of a certain universality and robustness of these 
features. 

9.1.1. Mass gap Since the critical solution is one of finite mass M c > M g > we are 
dealing with a type I phase transition in phase space. This is in contrast to the type 
II critical phenomena in GR |27] . It is interesting that the collapse of a massive scalar 
field in spherical symmetry does also exhibit type I phenomena |12]. However, unlike 
here, there continue to exist type II phase transitions for a sufficiently small scalar field 
mass. The fact that for a massive scalar field there are both type I and type II critical 
phenomena further strengthens the argument we presented in Section HI according to 
which no scalar field potential can be added to the matter Lagrangian to reproduce the 
modified field equations we are solving. 

9.1.2. Oscillations We would like to emphasize that the oscillations are not a numerical 
artefact. This is because they show up for both initial scalar field profiles in (1321) and 
( I33jl . as well as for all the functions q that we implemented (see Section f9. 1.5ft . In the 
case of GR, the presence of oscillations in the mass scaling law is well known |41[ I4TJ] 
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Equation 




©2D 


©5} 


Domain 


[-14,-6] 


[-14,-6] 


[-12,-4] 


[log 10 (M A H)]min 


-0.160 


-0.162 


-0.162 


5[log lo (0 o - <f>o)] 


1.342 


1.342 


1.342 


5[log 10 (M AH )] 


0.216 


0.219 


0.217 



Table 1. A comparison of the specifics of the modified collapse for different initial 
data (|33|) and a variation of w (see (|35|) ) instead of <po (see (|32p). 

and its period in log 10 -log 10 plots is approximated by 1.998. The period of the GR 
data in Figure [3] is about 2.010 (we used an inverse discrete Fourier transformation over 
the interval [—13,-3] on the horizontal axis). Further comments on the oscillations 
exhibited by the modified data are made in the concluding Section [TQl 

9.1.3. Overlap The fact that the points from the modified collapse agree with those 
from the classical collapse for a range of "large" amplitudes was to be expected. It 
confirms the intuition that energetic initial data configurations behave as if the spacetime 
is classical; the effective quantum corrections cannot be "seen". Note that in Figure [3] 
there are no points to the right of log lo (0o — 4>q) — as this is where the initial data is 
so "heavy" that it basically contains an apparent horizon. 

9.1.4- Scaling with A Since Mah is dimensionless we can read off Figure [3] the impact 
a change in A has. If we decrease it, the mass Mah = AMah of the apparent horizon 
has to become smaller and the data points corresponding to Mah & re translated down 
towards larger negative numbers. However, the GR points serve as fixed points in this 
translation ("attractors") in the sense that the overlapping region becomes larger. For 
A — > the interval over which the modified data displays a mass gap disappears. That 
is, the oscillations are "stretched out" and the data is in accordance with the classical 
mass scaling law. 

9.1.5. Universality and robustness In this section we provide evidence for the claim 
that the mass gap, the oscillations, the overlap, and the scaling with A are 

1) universal with respect to the initial scalar field profile and 

2) robust under a change of the deformation function q. 

Universality As for 1), the plot of the classical and the modified mass scaling relation 
for the initial data specified in (133]) displays a behavior that is identical to that 
apparent from Figure [3] |35| . To be more precise, Table [TJ shows that the size of 
the mass gap (we measured it by the minimum value of log 10 (MAH) ; that is we 
determined [log 10 (M A H)]min) and the specifics of the oscillations (comprising the period 
<5[log lo (0o — <f>o)], which we determined by an inverse discrete Fourier transformation, 
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(a) A family of functions q. 



x 

< 



cuO 
O 



1 


-1 
-2 
-3 
-4 



qo = 9/27 



vw\ 

qo = 5/27 



GR 




-15 -13 -11 -9 -7 -5 -3 -1 1 

logioOo - 0o) 
(b) Mass scaling plots for various qo. 



Figure 4. A family of functions q parametrized by qo and defined in (|36|) is illustrated 
in Figure 4(a) Note that qo 



has to be less than or equal to 1/3 as otherwise the 
modified evolution equations are singular. The mass scaling plots are given in Figure 
|4(b)| As compared to qo = 9/27 and the GR case, there are five times more data points 
corresponding to qo = 5/27 in order to better resolve the minima of the oscillations. 



and the amplitude 5[log 10 (MAH)]) agree for both (132]) and f )33|) over the given intervals. 
To extend this notion of universality, we varied w and fixed 0o for the initial data ( 1321) . 
just like Choptuik did. Once again the features mentioned previously are present [35] 
and the specifics are the same. The details are given in Table [TJ 

Robustness Regarding 2), we considered three different one-parameter families of 



functions q. Two of them are discussed in [35J and Figure 4(a) shows the function 



q + (1 -3q )r + 3q r 2 -q r 3 , r < 1, 

q = { . , ( 36 ) 

r, r>l, 
for two values of the parameter 

qo = q(0). 

For this function we need to impose the restriction q < 1/3 to guarantee that the 
modified evolution equations are not singular due to a local minimum of q (see ( |27|) ). 
The function q in (136|) is interesting in that it reduces to q = r in the limit qo — > 0. 
Another aspect of ( )36l) is that it causes the characteristic directions c± to be spacelike 
for r < 1 and null otherwise. For r > 1 we can expect the causal structure of the 
spacetime to be identical to that implied by GR. 



Figure 4(b) shows that for two values of the qo in ( 1361) the mass gap, the oscillations, 
the overlap, and the scaling with respect to A are present features. The conclusion is 
that Figure H] gives evidence in favour of the robustness proclaimed in 2). 
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10. Discussion 

In this paper we derived a set of effective equations for the massless scalar field in 
spherical symmetry. These equations deform GR but are nonetheless such that the 
constraint algebra closes. The deformation function q can be interpreted as providing 
effective quantum gravity corrections based on an analogy with the inverse triad operator 
in LQG. Although this analogy is not exact, it is nevertheless of sufficient interest to 
numerically determine the resulting quantum effects on the gravitational collapse. 

An interesting fact leading to the final form of our Equations (126|) and (127|) is that 
the characteristics of the scalar field do not coincide with the null directions of the 
metric. In the deformed theory, the characteristics can be timelike, null, or spacelike, 
depending on the form of the corrections. They are null if and only if there are no 
quantum corrections. This is the central reason why the equations in double null metric 
coordinates cannot be put into a useful form. However, when coordinates adapted to 
the scalar field characteristics are used, the equations simplify dramatically. 

A physical interpretation of the mismatch between the characteristics and the null 
directions of the metric is that the deformation can introduce violations of the dominant 
energy condition depending on the choice of the functions Qu\ and thus q. This can be 
seen directly by noting that we can define (in the notation of (Jl|)) an effective energy 
density 

= 4«2< 2 )«f + QmHf) = (^) 2 «f + (*)'«?>■ 

Then, since the diffeomorphism constraint is not modified, the dominant energy 
condition becomes 

e> Mi_ 

A 

It appears that this inequality can be violated for a suitably chosen q. We can speculate 
that such a feature is common to all effective quantum gravity equations, however 
derived, because singularity avoidance without a violation of the dominant energy 
condition appears to be an unlikely occurrence. Our example provides one class of 
deformed and anomaly-free equations with this property. 

A violation of the dominant energy condition was also found by Ziprick and 
Kunstatter in [18J. They derive modified constraints from the dilaton action formulation 
of "our" model, because of which the corresponding Dirac algebra is necessarily anomaly- 
free. The question of how their approach relates to ours is an interesting one. 

In this paper we also numerically solved and analysed the modified collapse 
equations that we derived in Section HI The specific choice of q was guided by the 
conditions that its impact vanishes in the asymptotic region and that it regularizes the 
singular factors 1/r 2 in the gauge fixed Hamiltonian constraint. This is in agreement 
with the expectation that quantum gravity effects only modify regions of large energy 
densities. The same conditions for effective quantum gravity modifications for the 
gravitational collapse equations were also relied on in [T5l [T7] . 
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With the function q defined, we used a new combination of numerical methods to 
get an accurate code without any mesh refinement. That is, unlike in previous work 
[T3l [15], we use a uniform grid in all our simulations. What is special about this grid is 
that it discretizes the dynamical fields not along the null lines of the metric but along the 
characteristic lines of the evolution equations. In GR they of course agree but because 
of the modifications the characteristics are better suited in our case. 

The results of our analysis of the modified mass amplitude data are as follows. 
We found a mass gap, jagged oscillations, agreement with GR for sufficiently large 
amplitudes O > scaling with respect to the quantum length scale A, a universality 
pertaining to different initial data profiles, and a robustness under a change of the 
function q. What is interesting are the form of the oscillations and the robustness. The 
latter implies that the nature of the characteristic foliation defined by a particular q has 
no impact on the presence of the key findings but only on their fine structure. 

As for the oscillations, they are new in comparison to earlier work [15j on quantum 
gravity corrected collapse in double null coordinates. The origin of the oscillations must 
be either that we use not exactly double null coordinates or that the way we implemented 
the modifications q does not lead to a violation of the Dirac algebra. 

It is interesting that Ziprick and Kunstatter report in |43j (see also |44] by 
Taves and Kunstatter) of similar oscillations resulting from the classical collapse in 
Painleve-Gullstrand coordinates. In jTT] they show that their quantum gravity modified 
model leads to the oscillations too. However, their analysis is based on the very 
different Painleve-Gullstrand coordinates, which makes a comparison quite difficult. The 
Painleve-Gullstrand coordinates allow for an evolution past the formation of apparent 
horizons, whereas our code has to terminate when such a horizon forms. We think that 
while the appearance of the oscillations in Ziprick and Kunstatter's work is in some way 
a result of the Painleve-Gullstrand coordinates, the oscillations we find are more the 
result of the effective quantum gravity corrections q and the fact that they do not give 
rise to an anomaly in the Dirac algebra. 

A future project that is suggested by these considerations is a numerical analysis 
of our modified model in Painleve-Gullstrand coordinates. This should clarify what the 
effect of the closure of the constraint algebra is. Furthermore, in Painleve-Gullstrand 
coordinates it is also possible to determine the causal structure of the spacetime for 
various functions q. In this context, an interesting question is to which extent our 
model can be compared to that of a collapsing fc-essence [45J. Namely, it can be shown 
that the scalar field equation given in (123|) is equivalent to the equation □</> = 0, where 
□ = g^VuVu is defined by the effective metric 

ds 2 = -2d v rFdudv + q W. 

On the one hand, this metric g^ u implies that the scalar field evolves on a geometry, 
in which the smallest sphere has a surface area that is bounded from below by 
47r(min r > {q(r)}) 2 > 0. On the other hand, it defines a second causal structure that 
can be contrasted to the "background" structure given by g^. The latter will result in 



Modified general relativity as a model for quantum gravitational collapse 



20 



pictures that describe a fc-essence. Another interesting question is what is happening 
to the self-similarity when the modifications are in action. To answer it, we can refine 
our code by an implementation of Garfinkle's algorithm [13]. Finally, a related question 
that is still open is whether the critical solution of the modified collapse is an unstable 
star. 

Acknowledgments 

This work was supported in part by the Natural Science and Engineering Research 
Council of Canada. We thank Martin Bojowald and Juan Reyes for providing a 
double null form of the scalar field-gravity equations in connection-triad variables, 
which motivated this work. We also thank Jack Gegenberg and Gabor Kunstatter 
for discussions, and the Atlantic Computational Excellence Network for computing 
resources. 

References 

[1] R. Amowitt, S. Deser, and C. Misner. The Dynamics of General Relativity. In L. Witten, editor. 

Gravitation: An Introduction to Current Research. Whiley, New York, 1962. 
[2] J. A. Wheeler. Superspace and the Nature of Quantum Geometrodynamics. In C. M. DeWitt 

and J. A. Wheeler, editors, Battelles Rencontres: 1967 Lectures in Mathematics and Physics. 

Benjamin, New York, 1968. 
[3] B. S. DeWitt. Quantum theory of gravity. I. The canonical theory. Phys. Rev., 160(5):1113-1148, 

1967. 

[4] C. J. Isham. Conceptual and geometrical problems in quantum gravity. Lecture Notes in 
Physics, 396:123-229, 1991. K. V. Kuchaf. Time and Interpretations of Quantum Gravity. In 
G. Kunstatter, D. Vincent, and J. Williams, editors, Proceedings of the 4 th Canadian Conference 
on General Relativity and Relativistic Astrophysics, pages 211-314. World Scientific, Singapore, 
1992. C. Isham. Canonical Quantum Gravity and the Problem of Time. In L. A. Ibort and 
M. A. Rodriguez, editors, Integrable Systems, Quantum Groups and Quantum Field Theory, 
pages 157-287. Kluwer Academic, Dordrecht, Boston, 1993. 

[5] C. W. Misner. Minisuperspace. In J. Klauder, editor, Magic without Magic: John Archibald 
Wheeler, pages 441-473. W. H. Freeman, San Francisco, 1972. 

[6] M. Bojowlad. Loop quantum cosmology. Living Rev. Relativity, 11(4), 2008. 
|http : / / www . livingreviews . org/lrr- 2008- 4[ A. Ashtekar. Loop quantum cosmology: 
an overview. Gen. Rel. Grav., 41(4):707-741, 2009. 

[7] M. W. Choptuik. A Study of Numerical Techniques for Radiative Problems in General Relativity. 
PhD thesis, The University of British Columbia, 1986. 

[8] D. Christodoulou. Violation of cosmic censorship in the gravitational collapse of a dust cloud. 
Commun. Math. Phys., 93:171-195, 1984. 

[9] D. S. Goldwirth and T. Piran. Gravitational collapse of massless scalar field and cosmic censorship. 
Phys. Rev. D, 36(12):3575-3581, 1987. 
[10] D. Christodoulou. The problem of a self-gravitating scalar field. Commun. Math. Phys., 105:337- 
361, 1986. 

[11] M. W. Choptuik, D. S. Goldwirth, and T. Piran. A direct comparison of two codes in numerical 

relativity. Class. Quantum Grav., 9(3):721-750, 1992. 
[12] M. W. Choptuik. Universality and scaling in gravitational collapse of a massless scalar field. Phys. 

Rev. Lett, 70(1):9-12, 1993. 



Modified general relativity as a model for quantum gravitational collapse 



21 



[13] D. Garfinkle. Choptuik scaling in null coordinates. Phys. Rev. D, 51(10):5558-5561, 1994. 

[14] V. Husain and M. Olivier. Scalar field collapse in three-dimensional AdS spacetime. Class. 

Quantum Grav., 18(2):L1, 2001. M. Birukou, V. Husain. G. Kunstatter, E. Vaz, and M. Olivier. 

Spherically symmetric scalar field collapse in any dimension. Phys. Rev. Lett., 65(10):104036, 

2002. V. Husain. G. Kunstatter. B. Preston, and M. Birukou. Anti-de Sitter gravitational 

collapse. Class. Quantum Grav., 20(4):L23, 2003. 
[15] V. Husain. Critical behavior in gravitational collapse. Adv. Sci. Lett., 2(2):214-220, 2009. 
[16] V. Husain and O. Winkler. Flat slice Hamiltonian formalism for dynamical black holes. Phys. 

Rev. D, 71(10):104001, 2005. 
[17] J. Ziprick and G. Kunstatter. Dynamical singularity resolution in spherically symmetric black 

hole formation. Phys. Rev. D, 80(2):024032, 2009. 
[18] J. Ziprick and G. Kunstatter. Quantum corrected spherical collapse: a phenomenological 

framework. arXiv:1004.0525vl, 2010. 
[19] M. Bojowald, M. Kagan, P. Singh, H. H. Hernandez, and A. Skirzewski. Hamiltonian cosmological 

perturbation theory with loop quantum gravity corrections. Phys. Rev. D, 74(12):123512, 

2006. M. Bojowald and G. M. Hossain. Loop quantum gravity corrections to gravitational 

wave dispersion. Class. Quantum Grav., 24(18):4801, 2007. M. Bojowald and G. M. Hossain. 

Cosmological vector modes and quantum gravity effects. Phys. Rev. D, 77(2):023508, 2008. 
[20] J. D. Reyes. Spherically Symmetric Loop Quantum Gravity: Connection to Two- Dimensional 

Models and Applications to Gravitational Collapse. PhD thesis, The Pennsylvania State 

University, 2009. 

[21] B. K. Berger, D. M. Chitre, V. E. Moncrief, and Y. Nutku. Hamiltonian formulation of spherically 

symmetric gravitational fields. Phys. Rev. D, 5(10):2467-2470, 1972. W. G. Unruh. Notes on 

black-hole evaporation. Phys. Rev. D, 14(4):870-892, 1976. K. V. Kuchar. Geometrodynamics 

of Schwarzschild black holes. Phys. Rev. D, 50(6):3961-3981, 1994. 
[22] A. Ashtekar. New Hamiltonian formulation of general relativity. Phys. Rev. D, 36(6):1587-1602, 

1987. J. F. Barbero. Real Ashtekar variables for Lorentzian signature spacetimes. Phys. Rev. 

D, 51(10):5507-5510, 1995. G. Immirzi. Real and complex connections for canonical gravity. 

Class. Quantum Grav., 14(10):L144, 1997. 
[23] T. Thiemann. Quantum spin dynamics. Class. Quantum Grav., 15(4) :839, 1998. 
[24] M. Bojowald, T. Harada, and R. Tibrewala. Lemaitre-Tolman-Bondi collapse from the perspective 

of loop quantum gravity. Phys. Rev. D, 78(6):064057, 2008. M. Bojowald, J. D. Reyes, and 

R. Tibrewala. Non-marginal LTB-like models with inverse triad corrections from loop quantum 

gravity. Phys. Rev. D, 80(8):084002, 2009. 
[25] V. Husain and O. Winkler. Semiclassical states for quantum cosmology. Phys. Rev. D, 

75(2):024014, 2007. 

[26] P. A. M. Dirac. Lectures on Quantum Mechanics. Belfer Graduate School of Science, New York, 

1964. A. Hanson, T. Regge, and C. Teitelboim. Constrained Hamiltonian Systems. Accademia 

Nazionale dei Lincei, Roma, 1976. M. Henneaux and C. Teitelboim. Quantization of Gauge 

Systems. Princeton University Press, Princeton, New Jersey, 1992. 
[27] J. M. Martin-Garcia and C. Gundlach. Critical phenomena in gravitational collapse. Living Rev. 

Relativity, 10(5), 2007. http://www.livingreviews.org/lrr-2007-5 
[28] M. Berger and J. Oliger. Adaptive mesh refinement of hyperbolic partial differential equations. 

J. Comput. Phys., 53(3):484-512, 1984. 
[29] D. Christodoulou. Bounded variation solutions of the spherically symmetric Einstein-scalar field 

equations. Commun. Pure Appl. Math., XLVL1131-1220, 1993. 
[30] J. C. Strikwerda. Finite Difference Schemes and Partial Differential Equations. Society for 

Industrial and Applied Mathematics, Philadelphia, 2nd edition, 2004. 
[31] V. Husain and O. Winkler. Quantum resolution of black hole singularities. Class. Quantum Grav., 

22T127-L134, 2005. 

[32] H. Halvorson. Complementarity of representations in quantum mechanics. Stud. Hist. Philos. Sci. 



Modified general relativity as a model for quantum gravitational collapse 



22 



B Stud. Hist. Philos. Mod. Phys., 35(l):45-56, 2003. A. Ashtekar, S. Fairhurst, and J. L. Willis. 
Quantum gravity, shadow states, and quantum mechanics. Class. Quantum Grav., 20(6): 1031, 
2003. 

[33] J. D. Hoffman. Numerical Methods for Engineers and Scientists. CRC Press, New York, 2nd 
edition, 2001. 

[34] M. W. Choptuik. Consistency of finite-difference solutions of einstein's equations. Phys. Rev. D., 
44(10):3124-3135, 1991. F. Pretorius. Numerical Simulations of Gravitational Collapse. PhD 
thesis, The University of British Columbia, 2002. F. Pretorius and M. W. Choptuik. Adaptive 
mesh refinement for coupled elliptic-hyperbolic systems. J. Comput. Phys., 218(l):246-274, 
2006. 

[35] A. Kreienbuehl. Quantum Cosmology, Polymer Matter, and Modified Collapse. PhD thesis, The 

University of New Brunswick, 2011. 
[36] J. Thornburg. Adaptive mesh refinement for characteristic grids. Gen. Relat. Gravit., 43(5):1211- 

1251, 2011. 

[37] J. M. Bardeen and T. Piran. General relativistic axisymmetric rotating systems: coordinates and 
equations. Physics Reports (Review Section of Physics Letters), 96(4):205-250, 1983. 

[38] M. Alcubierre. Introduction to 3+1 Numerical Relativity. Oxford University Press, Oxford, 2008. 

[39] S. Ayal and T. Piran. Spherical collapse of a massless scalar field with semiclassical corrections. 
Phys. Rev. D, 56(8):4768-4774, 1997. 

[40] C. Gundlach. Understanding critical collapse of a scalar field. Phys. Rev. D, 55(2):695-713, 1997. 

[41] S. Hod and T. Piran. Fine structure of Choptuik's mass-scaling relation. Phys. Rev. D, 
55(2):R440-R442, 1997. 

[42] P. R. Brady, C. M. Chambers, and S. M. C. V. Goncalves. Phases of massive scalar field collapse. 

Phys. Rev. D, 56(10) :R6057-R6061, 1997. 
[43] J. Ziprick and G. Kunstatter. Numerical study of black-hole formation in Painleve-Gullstrand 

coordinates. Phys. Rev. D, 79(10):101503(R), 2009. 
[44] T. Taves and G. Kunstatter. Higher dimensional Choptuik scaling in Painleve Gullstrand 

coordinates, 1997. arXiv:1105.0878vl. 
[45] R. Akhoury, D. Garfinkle, and R. Saotome. Gravitational collapse of fc-essence, 2011. 

arXiv:1103.0290vl. D. Leonard, 2011. talk at the Black Holes VIII conference. 



